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FOREWORD 

The  objective  of  this  program  is  to  develop  general  procedures  for  the 
determination  of  electromagnetic  scattering  from  arbitrary  surfaces  with  appli¬ 
cation  to  optimization  of  vehicle,  antenna,  and  radome  designs. 

This  document  is  principally  a  compilation  of  the  following  technics  reports 
and  memos  issued  during  the  course  of  this  investigation. 

1.  Technical  Data;  ’’Scattering  of  a  Plane  Magnetic  Wave  from  a 
Conducting  Circular  Cylinder"  -  Fred  Oshiro,  1  February  1963. 

z.  Technical  Data:  "Scattering  of  a  Plane  Electric  Wave  from  a 
Conducting  Circular  Cylinder"  -  Fred  Oshiro,  7  February  1963. 

3.  Technical  Data:  "A  Source  Distribution  Technique  for  the  Solution 
of  General  Electromagnetic  Scattering  Problem"  -  February  1963. 

1.  NB63-125;  "Proposal  for  a  Source  Distribution  Technique  for  the 

Solution  of  General  Electromagnetic  Scattering  Problem"  -  April  1963. 

5.  Addendum  &  Technical  Data  (Item  3)  dated  23  May  1964. 

6.  Preliminary  Data,  "A  Source  Distribution  Technique  for  the  Solution  ot 
General  Electromagnetic  Scattering  Problems",  dated  14  May  I ’>64. 

7.  Addendum  to  Item  (6)  dated  14  December  1964. 


The  technique  described  herein  is  being  further  developed  and  extended  td 
solve  arbitrary  geometries  unde;*  contract  AF  33(615)-  jlb6  ••  "Calculation  of 
Radar  Cross  Section"  -  A  right  Patterson  Air  Force  Base,  Ohio 
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SUMMARY 

A  Source  Distribution  Technique  ("SDT")  for  the  solution  of  general  electromagnetic 
problems  has  been  developed  by  the  Electronic  Systems  Research  Group  of  Northrop 
Nor.s  ir.  his  technique  embodies  a  comparatively  simple  and  unrestricted  computa-  ~ 
tnmal  method  applicable  to  the  complex  problem  of  scattering  from  arbitrarily  shaped 
geometries.  The  procedure  exhibits  a  greater  degree  of  flexibility  and  accuracy  than 
present  day  approximation  techniques.  It  is  a  numerical  procedure  for  determination 
o!  tlie  radar  c  ,jS  section  for  complex  targets,  based  on  an  integral  equation  approach. 

It  cun  be  employed  to  determine  the  electromagnetic  scattering  from  arbitrarily  shaped 
two-and-three-dimensional  geometries  in  the  difficult  resonant  frequency  ranges  in 
which  classical  solutions  are  not  available.  Combined  with  the  Norair  Physical  Optic., 
procedure  the  ”SDT"  can  be  employed  to  determine  the  electromagnetic  scattering  of  a 
variety  of  shapes.  Including  large  scattering  objects;  t.e.  >2  A. 

Verification  of  the  technique  was  performed  on  the  Norair  radar  reflectivity  range. 
An  image-plane  was  utilized  to  measure  the  current  distribution  and  scattered  fields. 

A  brief  outline  of  the  basic  procedure  is  described  together  with  a  few  examples  to 
demonstrate  its  validity. 
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SYMBOL  DEFINITIONS 


a  =  arbitrary  constant  vector 

A 

E(r)  =  electric  fie  ld  vector 

F  =  unknown  current 

$(r)  =  induced  current  vector 

Ft  =  C  +  JD  =  transverse  current 

F  =  A  +  JB  =  longitudinal  current 

Z 

G(r',r  )  =  Green's  function 

ft(r)  =  magnetic  field  vector 

H0(n(kr)  *  aero  order  Henkel  taction  of  file  first  kind 

H|^(kr)  »  first  order  Hankel  taction  of  the  first  kind 
J  -  tat  imaginary  numbs r  (J2  «  -1) 

k  *  ^  »  wave  anther 

K  «  incident  field 

n  =  unit  normal  vt  .tor  -Jtx  +  my  +  ni 

np  -  unit  normal  at  point  F  *  + rapy  + tip  z 

n  =  number  of  surface  elements 

2  2  2 

r  -  k  +  y  +  t  =  observation  point  in  region  V 

S  =  surfaces  of  region  V 

t  -  transverse  component 


NOTATIONS 


Subscripts: 

P  observation  point  on  acatterer 


Q  _  source  point  on  scatterer 

1  -  obstacle  surlace 

2  -  outer  surCace 

t  =  transverse  component 
l  =  longitudinal  component  in  z  direction 


Superscripts: 

primed  variables  *  represents  source  points 
unprimed  variables  *  represents  observation  points 
bu.  « incident  wave 

sc 


-  scattered  wave 


1.0  INTRODUCTION 
1. 1  Problem  Area 

Rapidly  advancing  weapon  system  technologies  have  created  a  renewed  interest 
in  the  study  of  scattering  and  diffraction  probleo-s.  A  current  concern  is  with 
evaluation  and  control  of  radar  cro.ss  section  (RCS).  Much  effort  has  been  directed 
at  reducing  RCS  by  camouflaging  techniques  aid  by  use  of  radar  absorbing  material 
(RAM);  additional  efforts  are  directed  at  RCS  enhancement  of  drones  and  decoy 
applications  for  penetration  aids.  There  is  an  acute  interest  in  reflectivity  signature 
for  vehicle  identification.  Northrop  Norair’s  studies  in  this  field  have  beeu  motivated 
by  a  need  to  op*  nise  vehicle  shapes  compatible  with  both  aerodynamic  performance 
and  RCS  requirements.  Solutions  to  these  scattering  problems  are  of  vital  interest 
to  the  design  engineer. 

The  computation  of  scattering  from  an  arbitrary  obstacle  is  intrinsically  compli¬ 
cated,  since  it  depends  mi  several  parameters  as  frequency,  polarisation,  direction 
of  incident  wave,  observation  angle,  configuration,  and  material  of  the  body.  The 
fundamental  problem  in  computation  Is  die  explicit  determination  of  the  far  field 
parameters  of  amplitude,  phase,  and  polarisation  by  solution  of  the  basic  wave  equation. 

TracKtionally,  solutions  to  the  wavs  equation  have  been  achieved  by  the  separation 
of  variables  or  eigenfunction  expansion,  asymptotic  evaluation  of  an  integral  equation, 
and  the  optical  approximation  techniques.  The  restrictions  imposed  by  these  methods 
have  limited  their  usefulness  in  many  practical  engineering  problems.  The  require¬ 
ment  that  the  body  configuration  be  completely  defined  in  a  coordinate  system  has  led 
to  special  mathematical  functions  that  are  tabulated  for  only  a  few  simple  shapes. 
Asymptotic  expansions  lead  to  infinite  series,  whose  convergence  behaviors  are  deter¬ 
mined  by  the  "electromagnetic  else"  of  the  body,  therefore,  for  practical  reasons,  only 


"small"  shapes  have  been  examined.  In  contrast,  the  "optical"  techniques  are  valid 
only  for  'large"  shapes. 

An  immediate  dilemma  arises  when  the  design  engineer  is  faced  with  a  scatter¬ 
ing  problem  in  which  the  geometrical  shape  does  not  "fit"  the  above  requirements.  The 
problem  usually  is  approached  by  intelligent  application  of  known  solutions  from 
simpler  shapes,  or  by  relying  heavily  on  reflectivity  range  data.  A  need  then  exists 
for  a  general  method  for  evaluating  the  scattering  and  diffraction  problems  encountered 
in  the  design  of  practical  electromagnetic  systems. 

1.2  state  of  the  Art 

Since  the  solution  of  electromagnetic  scattering  tnd  diffraction  problems  is 
difficult,  the  introduction  of  the  above  restrictions  was  not  without  foundation.  In 
order  to  gain  insight  as  well  as  attach  physical  significance  to  the  theoretical  studies, 
it  was  natural  that  the  classical  atuuies  be  concentrated  on  small,  simple,  and  idealized 
bodies.  Most  of  the  present  day  approximation  procedures  for  complex  structures  are 
usually  tested  against  the  exact  solutions  from  simpler  shapes.  The  simple  body 
solutions  provide  the  only  means  for  calibrating  experimental  techniques  to  measure 
scattering  and  diffraction  phenomena.  Also  the  knowledge  gained  by  a  detailed  study 
of  the  simpler  shapes  is  used  as  stepping  stones  to  the  study  of  more  complicated 
shapes. 

The  accomplishments  of  Siegel  and  his  group  at  the  University  of  Michigan 
together  with  the  Antenna  Laboratory  Group  at  Ohio  State  University,  constitute  a 
major  contribution  in  recent  years  (References  1  and  2).  Their  procedure  was  to 
assume  that  the  various  parts  of  a  large  complex  shape  (such  as  an  aircraft)  scatter 


independently.  They  approximated  the  parts  by  relatively  simple  shapes,  computed 
the  scattering  from  each  shape,  and  combined  the  individual  contributions  to  obtain 
the  total  scattered  iield.  Their  results  yield  good  agreement  between  theory  and 
experiment.  However,  for  critical  aspect  angles  large  errors  arise  which  are 
generally  attributed  to  the  Inability  of  the  method  to  properly  model  the  resonance 
of  the  fine  structure  (i.e. ,  ducts  and  engine  nacelle). 

In  attempting  solutions  to  more  general  configurations,  the  emphasis  is  usually 
placed  on  the  desirability  of  attaining  purely  analytical  solutions.  The  study  of 
numerical  solutions  has  been  neglected  but  not  overlooked,  hi  1959  Sinclair  (Refer¬ 
ence  3)  emphasizes  the  need  for  more  general  methods  of  solving  boundary-value 
problems.  He  also  pointed  out  that  there  may  be  certain  advantages  of  numerical 

solutions  in  the  "resonance  region. "  The  usefulness  of  numerical  techniques  in  the 
solution  of  potential  flow  problems  (Laplace's  equation)  in  aerodynamics  and  hydrody¬ 
namics  has  been  demonstrated  by  Smith  (Reference  4)  and  Hess  (Reference  5).  In 
addition  to  solving  general  problems  involving  two-dimensional  and  axially  symmetric 
bodies,  they  have  provided  general  techniques  for  the  solution  of  arbitrary  three- 
dimensional  shapes.  There  is  a  very  limited  supply  of  numerical  procedures  for 
general  electromagnetic  scattering  problems. 

It  is  the  purpose  of  this  report  to  describe  a  method  based  on  the  numerical 
solution  of  an  integral  equation  which  is  capable  of  giving  excellent  results  in  this 
range.  Additionally,  by  application  of  the  principles  of  physical  optics,  the  technique 
is  extended  to  include  "large"  objects  and  has  the  nique  feature  of  eliminating  much 
of  the  guesswork  associated  with  the  current  customary  physical  optics  procedures. 
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In  the  volume  V  enclosed  by  surfaces,  Sj  and  assume  there  exists  two 
scalar  potential  functions  ^  (r)  and  Q  (r',  r)  such  that  #(r)  satisfies  the  homo¬ 
geneous  wave  equation, 


(V2  +  kJ)  *  (r)  •  0 

^(r)-»0as  r-«« 


(l) 


and  the  Green's  function  G  (r\  r)  satisfies  the  wave  equation. 

(V2  +  k2)  G  (r»,  r)  «  -  6  (r  -  r’), 
G  (r*,  r)-0asr— 


4  - 


00  ♦ 


(2) 


Applying  Green'*  Theorem  to  $(v)  and  G  (r',  r)  we  obtain 


/  |c  (r'(  r)  IP’2  +  k2)  1/1  (r ')-  0  (r'l  (v'2~  k2)  0  (r'  ,  rjjdv' 
/  jofr’.rtjir  (Mr')- *(r')jSr  0(r’,r)|<l»' 

S,  +  Sj  • 


where  the  normal  derivative  on  die  aurface  ia  taken  in  the  direction  of  the  inward 
normal. 

However ,  in  virtue  of  Equation  1, 

/  |  G  <r',  r)  (r*2  +  k2)  *  (r*)|  d»'  -  0  ,  (4) 

also  from  the  definition  of  the  delta  function 


-  /  |  *  <r')  (f*  ♦  k1)  O  <r',  r))dv». 

(5) 

-  /  #  <**•)  |  -  4  ( r  -  r*)  ]  dv»  -  *  (r)  . 


s  * 


(6) 


The  total  field  at  r,  is  composed  of  an  incident  and  scattered  field. 


*toU1<r).*‘nc(r)  +  **c<r). 


&  can  be  shown,  (References  6  and  7)  that  the  integral  over  Sg  reduces  to 
iuc 

< U  (r)  and  the  integral  over  represents  the  scattered  field. 

Equation  6  becomes 

*toM  <r>  =  *lnc  (r)  +  /  |  O  <r\  r)^-  *  <r’)  -  *  <r')  jgr  G  (r1,  r)  |  ds' 

S1 

(8) 

and  represents  the  integral  equation  of  the  field  at  r  in  terms  of  an  incident  field 
contribution  and  a  surface  integral  contribution. 

2.1.1  l  IXMssk  (E..100)  ProcmUoa 

For  the  case  of  TM  mode  propagation  in  the  two  dimemioMil  domain  the  electric 
field  can  be  expressed  as 

I(r)- #(#,•)»-*,  (t)i  (9) 

where  s  is  taken  as  the  separation  direction  and  t  represents  the  transverse 
components.  (This  propagation  mode  corresponds  to  horixontal  polarization).  Sub¬ 
stituting  this  expression  in  Equation  8  yields, 


txt<m  (t)  -  Etlnc  (t)+  /  |g  (t\  t)j£r  (t>) 

81 


-  Eztotal  (t')^r  G  (f,  t)  |  d»' 


*  8  - 


The  free  space  Green's  function  In  the  two  dimensional  domain  is  the  sero 


order  Hankel  function  of  the  first  kind, 

G(P*,  P,  H0(1)(kVp2+^’2-  2PP»cos  (d-d')  (U) 

<*  G(t»,  t)*^  HQ(1)(k|t  -  t'|)  ^  (12' 

If  Sj  is  perfectly  conducting,  the  total  electric  field  vanishes  on  the  surface, 

£  totai  (t)  M  E  inc  (t)  +  £  ac  (t)  =  0  (13; 

z  z  •  • 

Therefore  when  the  observation  point  is  taken  on  the  surface.  Equation  10 
reduces  to 

0  •  Eltoc  <t)  +  J  f  H0(1)  (k|t  -t’D-j^r  EztOU1  (»')  ds'_  (14) 

81 

The  quantify  (t*)  Is  the  unknown  current  induced  on  by  the 

incident  electric  field.  The  direction  of  current  flow  is  in  the  direction  of  the  electric 
field,  e.g. ,  the  s  direction.  Setting, 

K  (t’>  ■  &  E,totil  (t,)  *  (15) 


and  introducing  subscripts  P  and  Q  for  the  observation  and  source  points  respectively 
on  the  surface  of  the  scatterer,  we  obtain  the  desired  integral  equation  for  the 


unknown  current  F„ 


(1q)| 


0  -  E2tac  (tp)  +  f-  /  [  H^1*  (k  |tp-  tQ| )  Fz  (tQ)  |  d,  Q  (18) 

* 
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2. 1.1.2  IE  Mode  <Hx  tnc)  Propagation 

For  the  case  of  TE  mode  propagation,  the  magnetic  field  vector  may  be 


expressed  as  follows: 


H  (r)  =  (t>,  B)  2  =  Hz  (t)  z 


where,  again,  z  is  taken  as  the  separation  direction.  (This  mode  of  pro¬ 
pagation  also  corresponds  to  vertical  polarization). 

Substituting  this  expression  In  Equation  8  yields, 

H^U>U1  ft)  „  H^inc  (t)  +  /  |Q  (t.(  Hztotal  (f) 

S1 

„  total  Ii|\  0  n  ft t  *\  1  rf.t 


-  Hst0U1  (t*)^  G  <t»,  t)  |  da* 


The  boundary  condition  on  Sj  for  die  TE  case  is 

?JL  (t>)  ■  o  (is) 

If  the  observation  point  is  taken  on  8^,  then  Maue  (Reference  8)  has  shown  that 
the  induced  current  at  the  point  of  incidence  is  the  average  value  of  the  magnetic  field 
ut  that  point,  consequently 


H  (t) 


=  Hzlnc  <t)  -  /  Hstotal  (t'>5*r  Ho(1>  Qc  It  -  t*| )  d»'#  (20) 


The  quantity  Hztotal(t')  is  the  equivalent  unknown  current  induced  on  Sj  by  the 


incident  magnetic  field  and  flows  transverse  to  the  z  direction. 


F.  (t*)  =  H  toUi  (t»)  t 

v  Z  • 
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The  Integral  equation  i«, 


H*t0U1  V  inc 
S—i - —  "  H,  (tp) 


<kitP*tQl)‘  Ft  V  c°s("> 


(2?; 


2. 1.2  scattering  Cross-Section 

The  scattered  fields  may  be  evaluated  from  the  current  distribution  by  applica¬ 
tion  of  Kirchoffs  integral 

tfSC  (r)  «  -  /  p*  G  (r\  r)xF  (rf)  ds*  (23) 

S 


Since  the  observation  distance  is  usually  taken  to  be  very  large,  the  asymptotic- 
forms  of  the  Hankel  functions  may  be  introduced  for  the  Green's  function  to  give: . 

1/2 

S,"  (r)  -  (^)  sJ  0tr  -  n/i)  / Ft«lkr'  c“  <•  ’  •’>  cos  (ft,  *-)  ds-  (24) 


Ef“  (r)  - 


r/4) /F.eJkr-cos(«-  s-)d>1  @5) 

s 


The  radar  cross  section  of  two-dimensional  cylinders  is  defined  as  the  power 
scattered  per  unit  axial  length  to  the  incident  power  and  has  dimensions  of  length. 
Thus, 


where  the  expressions  E.*6^)  and  H  “(r)  are  the  scattered  fields  from  Equations 

Z  Z 

24  and  25  respectively  and  foe  Integration  is  taken  over  the  closed  contour  of  the 
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scattcrer.  For  simplicity  the  scattering  cross-section  has  been  normalized  to 
unit)-  incident  power. 


2. 1.3 


ImerTcH 


The  use  of  an  integral  equation  to  formulate  a  scattering  problem  is  well  known 
and  has  been  the  basis  of  many  approximation  techniques.  (References  8-14).  The 
present  technique  differs  from  classical  treatments  in  that  numerical  methods  are 
exploited  to  obtain  direct  solution  for  the  unknown  quantities  rather  than  formal  indirect 
methods.  The  solution  of  Equations  16  and  22  ,  for  example,  may  be  achieved  by 
first  approximating  the  contour  Sj  by  n  straight  line  segments  and  then  setting  the 
current  distribution  constant  over  each  interval. 

Fxl*At  +  1Bl  (27) 

Ftt  *  ci  +  J  Dl  (28) 

This  assumption  is  valid  only  if  the  line  segment  is  less  than  .  12X.  Inserting 
this  expression  in  Equations  16  and  22  respectively  yields 


v  - 


■H2inC  (tp)  = 


|Cj(tQ)  +  )Pt(tQ) 


+  V"  I  S.  +  ^  Di  J  *P  ”  *Q  P  cos  (n*  "  *q))  ^8< 


The  integration  of  the  Hankcl  function  is  carried  out  numerically  via  Simpson's 
rule.  The  integral  equation  is  now  recast  into  a  set  of  finite  difference  equations, 
wherein  the  complex  coefficients  of  the  unknown  currents  form  an  n  x  n  matrix,  n 
being  the  number  of  sub-surfaces  describing  the  body.  The  excitation  function  is  the 
impinging  wave.  The  unknown  currents  are  determined  by  solving  a  matrix  equation 
of  the  type 

[t]  H  -  [«]  (3.) 


where 


-  unknown  currents 

=  Kernel  of  the  integral  equation 

-  incident  field  , 


Applications  of  matrix  inversion  techniques  lead  to  a  direct  solution  of  the 
currents. 


2.1.4  Application  to  Two- Dlmori tonal  Geometries 

To  demonstrate  the  validity  of  this  technique  a  classical  problem  is  examined, 
l.e. ,  scattering  from  an  infinite  circular  cylinder.  The  problem  is  specialized  for 
normal  incidence  of  either  a  vertically  or  horizontally  polarized  wave,  resulting 
in  either  transverse  or  long  tudinal  type  currents  respectively.  These  simplifications 
were  introduced  to  provide  suitable  comparisons  with  available  classical  results. 
Figures  1  and  2  are  comparisons  of  the  current  amplitude  and  phase  with  the  classical 
series  expansion  results.  (Inference  15).  In  both  examples  the  diameter  of  the 
cylinder  was  chosen  to  be  X/2  ;  and  twelve  points  were  distributed  equidistant  around 
the  periphery  tr  carry  t  the  numerical  integration. 
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The  bistatic  scattering  cross-sections  of  the  circular  cylinder  for  vertical 
and  horizontal  polarization  arc  shown  in  Figures  3  and  4  respectively.  The  trans¬ 
mitter  is  fixed  at  the  180  degree  position  and  the  radar  cross  section  is  observed  at 
various  histatic  angles  The  quantity  characteristics  of  energy  scattered  back  toward 
the  source  is  commonly  known  as  the  monostatic  cross  section  of  the  obstacle  where¬ 
in  the  transmitter  and  receiver  are  both  located  at  the  same  point.  The  forward  or 
total  scattering  cross-section  is  the  value  at  0  degrees.  Where  possible  theoretical 
results  have  been  introduced  for  suitable  comparison  with  the  "SDT"  data.  Normali¬ 
zation  of  the  results  is  accomplished  by  multiplying  the  cross-s action  by  the  wave 
number,  k.  (The  designation  ‘’SDT”  in  quotation  is  used  to  distinguish  this  technique 
from  other  proce-’res  developed  during  the  course  of  this  study). 

The  flexibility  of  the  "SDT"  for  the  solution  of  non-circular  geometries  w  ill  be 
demonstrated.  Earlier  contributions  to  generalize  the  circular  cylindrical  scattering 
problems  to  non-circular  shapes  were  made  by  Fock  (Reference  16).  His  "universal 
function"  enabled  him  to  determine  the  surface  currents  in  the  penumbra  region  near 
to  and  Including  the  shadow  boundary.  More  recently,  Wetze!  (Reference  17)  follow¬ 
ing  the  lead  of  Fock  derived  a  new  set  of  functions  to  obtain  the  currents  in  the  deep 
shadow  region  where  Fock's  approximation  was  not  adequate.  Wetzel  applied  his 
functions  for  the  solution  of  the  current  distribution  on  a  perfectly  conducting  elliptic 
cylinder  and  confirmed  his  theory  with  measured  data.  The  "SDT"  was  applied  to 
the  same  geometry  for  comparison.  As  shown  in  Figures  5  and  6  the  "SDT"  shows 
still  greater  correlation  with  the  measured  data,  especially  in  the  deep  shadow  region. 
Its  blstattc  scattering  cross-section  is  shown  in  Figure  7.  The  sharp  peak  in  the 
forward  direction  is  characteristic  of  large  scatterers. 

Next,  the  scattering  cross-section  of  two  different  airfoil  sections  shown  in 
Figures  8  and  9  were  considered.  An  incident  magnetic  plane  wave  propagating  along 
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the  x-axls  Impinges  on  the  surface  of  the  perfectly  conducting  two  dimensional 
geometries.  Classical  results  to  these  geometries  have  not  been  possible  because 
of  the  complexity  of  solutions.  The  airfoil  shape  was  approximated  by  49  straight 
line  segments  distributed  along  the  periphery  and  the  "SOT"  was  applied.  The  bi- 
etutic  cross-sections  are  shown  in  Figures  10  and  11. 

Several  important  conclusions  concerning  the  behavior  r  these  scatterers 
can  be  made  which  otherwise  would  be  difficult  if  not  impossible  to  predict  by 
approximation  techniques  or  by  intuition. 

1.  A  large  back  lobe  is  developed  at  150  degrees  in  the  direction  of  the 
trailing  edge. 

2.  Although  airfoil  section  No.  2  is  smaller  than  No.  1,  its  major  back  lobe 
is  larger. 

3.  The  monostatic  cross  sections  (at  180°)  of  the  airfoil  sections  are 
approximately  equal. 


Aside  for  the  additional  complexity  of  vector  fields,  the  derivation  of  the  inte¬ 
gral  equation  is  similar  to  the  two  dimensional  case. 


In  the  volume  V  enclosed  by  sur&ces  8^  and  S2  assume  there  exists  two  vector 
A  A  A 

fields  H(r)  and  G  (r',  r)  such  that  H  (r)  satisfies  the  homogeneous  vector  wave 


equation 

Vx  V  x  H  (r)  -  k2  H  (r)  =  0 
H  (r)— ■-Oas  r  — x 


(32) 
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>**  rfytawsw"  ’ 


...» 


H  (r)  =  Magnetic  Field 
G  (r;  r  )  =  <P  (r;  r  )  a 

,  e  i  k  lr"r'l 
“,r-  r>  =  4ff|r^r 


I  rp  “  rQ  ! 


V  Si 


A 

A  H(F)/ 
i  G(r\  r)  / 


and  the  vector  Green'*  function  G(r',r)  satisfies  the  Inhomogeneous  wave  equation 


V  x  Fx  G  (r’,  r)  -  k2  G  (r'»  r)  =  -  I  (r  -  r*)  a 


G  (r*,  r)— 0  m  r  — •-* 


Applying  Green's  theorem  to  the  wave  equation  we  find 

/  |g  (r',r)*^F’x  F’x  +  k2J  H  (r’)  -  H  (r1)’  |V'  x  F’  x  +  k2J  G  (r%  r)  J  d  v' 

.  A  -  (34) 

-J  [h  (r1)  x  r  x  G  (r'»  r)  -  G  (r*,  r)  x  F*  x  H  (rf)  I  •  n’ds' 


V*2 


where  the  normal  n*  is  directed  inward  toward  the  volume  V . 


The  left  hand  side  of  Equation  34  reduces  to 


-  /  H  (r')  •  [-  4  (r  -  r’)  a]  dv'  --  H  (r)  .  a  c  (35) 

The  vector  *a  is  assumed  constant  and  arbitrarily  oriented  within  volume  V. 
Equation  34  becomes, 

H  (r)  •  a  =  /  j  £n  x  P’x  H  (r’)  J  •  G  <r\  r)  +  nxH  (r*)  •  |V  x  G  (r\  r)  J  J  d  a 


■/if 


x  V'  x  H  (r')  1  •  G  (r’» r.i  +  nxH  (r*)  •  V'  x  G  (r*, r)  j  d  s'  (36) 


Substituting  H  (r)  ^(r)  -ii80  (r)  into  the  first  integral  of  Equation  36 


/  =  -/  j  n  x  [r*  x  H1"®  («•*)  j  •  G  (r’,r) 


S2  S2 


+  [nx  H111®  (r*)]  *  [r*  x  G  (r',  r)]  j  d  s'  (37j 
^  jn  x  [r  x  H1®  (r*)]  •  6  (r't  r)  +  [n  x  H1®  (r*)]  •[*•*&  (r’,  r)]  [  d  s’ 


Applying  Green's  theorem  in  reverse  to  the  first  integral 


-  /  Jnx  Tr  xHln®  (r’)j-  G(r’,r)+  £nxHinc(r’)]  •  [V  x  G  (r’,  r)]  j  d 
S2 

=  [  [6  (r1,  r)  •  (F*  x  7’  x  +  k2)  Hinc  (r’) 

V  (38) 

-  H1”®  (r’) .  (F‘  xf'x  +  k2)  G  (r’,  r)  J  d  v’ 
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/  H*110  (r’)  [-  6  (r  -  r’)  a]dv'  =  Hinc  (r)  • 


(3») 

Continued 


Ag«  A v 

The  second  integral  vanishes  identically  since  H  (r')  and  G  (r *,  r)  both  satisfy 
the  radiation  condition.  Therefore,  Equation  36  reduces  to 


H  (r)  •  a  =  Hlnc  (r)  •  a 

-  /  j^n  x  7’  x  H  (r‘)  J  •  G  (r',  r)  +  n  x  H  (r')  •  ^'xG  (r’,  r)J  |  d  s' 


And  from  Maxwell's  equation, 


V  x  H  (r')  ■  J  w  «  E  (r*^ 


we  find 


H  (r)  •  a  =  Hlnc  (r)  •  a  -/  j  J««  £nxl  (r') J  •  G  (r’, r) 


+  [n x  H  (r*)]  'pxO  (r*, r)]  J  d  s’ 


If  obstacle  S1  is  assumed  perfectly  conducting  then  the  tangential  component 

A  A 

of  the  electric  field  vanishes,  n  x  E  'r')  =  0  and  additionally  the  current  Induced  on 
Sj  Is  given  by  n  x  H  (r *>  =  F  (r'). 

Equation  41  may  be  rewritten  as 
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H  (rp)  •  a  =  Hinc  (rp)  •  a  -  /  F  (r^)  •  [VQx  G  (rQ  ,  Tp)J  d  sQ. 


x  G  (tqj  rp)=  FQx^(rQ,  rp)a  =  «p(FQxa) 

+  ^  ^rQ>  rp)  x  a  *  V  q  V  Uqi  fp)  X  a  since  a  is  a  constant  vector. 

Hence, 

H  (rp)  •  a  =  Hinc  (rp)  •  a  -  a  •  /  j  Fq  <P  (*q»  rp)  x  F  (r^)  J  d  sQ  (43) 

V 

If  the  observation  point  is  taken  on  the  surface  we  obtain  the  desired  integral 


equation 


F  (rp) 


—  *  +  Htac(rp)  x  np  -  *  <rQ>  rP)  x  F  Uq)  x  n  pJ  d  sQ 

Where  subscripts  P  end  Q  designate  points  on  the  surface  of  the  scatterer 

.  ,+  l*l  rp-  rQl 

*  (rQ>  V  *  T»|rp-PrQ|S 

9  (rp)  ■  current  at  point  P 

Vq0(tq  •  rp)  ■  gradient  of  0  with  respect  to  point  Q 


“lP-  * 

(rp-  fg)' 


p-  [1  +  Jk|rp-  rQ|]  <?(rQ,  rp)  »t;  ^  =  x,  y,  z 
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fHrp)  - 


-  incident  plane  wave  at  point  P 


np  -  unit  normal  at  point  P 

/A  A  A 

P  x  +  trip  y  +  np  z 

I  rp  ■  **qI  =  V(xp  '  xq>2  +  <yP  •  y q>2  +  -  zq)5 


S1  *  total  surface  of  the  scatterer 

dSj-j  =  element  surface  at  point  Q  . 
Vi 


2.2.2 


The  scattered  field  may  be  obtained  by  substituting  the  current  distribution  in 


Kirchoff's  integral. 


H8c(r)  = 


■it 


*  G  (r»,  r)  x  F  (r1)]  d  s 


r  r«+i*  *  I r  *  r’l  1  a  i 

■{[Vr - <£  +  Jk>*F<r’)Jd 


and  for  r  >  >  r* 


/[,!>“  lr  ‘  r'l*  F  (r')]d 

2 

The  scattering  cross-section  has  the  dimension  of  length  .  Thus 


•c  .  2  |H,C  (r)  |  2 

*  ’  4*r  IobTTTI 
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A 

where  H8C(r)  is  given  by  Equation  46  and  the  integration  is  carried  out  over  the  entire 
surface  of  the  scatterer. 


2.2.3  Numerical  Solution  of  the  Integral  Equation 

The  integral  Equation  44  represents  the  complete  formulation  of  the  boundary- 
value  problem  with  the  following  boundary  conditions: 

E  (r)  x  n  =  0  ,  H  (r)  x  np  =  -  F  (r) 


To  completely  evaluate  the  unknown  currents,  the  integral  equation  must  be  solved 
for  the  individual  components  of  $  (r).  For  the  spherical  geometry,  the  surface  may 
be  parameterized  in  terms  of  local  coordinates  9  and  o.  In  this  special  case  there  are 
two  components  of  the  unknown  current  and  the  governing  set  of  equations  are: 


Fg  =  2  sin  #  eJkaco#®+  8 ^-  /  d  if  d^^^-(l  +  jkr) 

0  -n  r 

j  £stn  9  sin  9  +  cos  9  cos  9  cos  -9)-  cos  -  #)J  F^ 


£  sin  ( ^  -  9)  (cos  9  -  cos  9)  J  F- 


2 


cos  9  cos 


-  1  ka  cos  6 

$  6 


(49) 


a _ 

“  F7T 


d  <P 


<i+)kr)S[]v[] 
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The  coupled  set  of  matrix  equations  are 


(50) 

(51) 


By  direct  extension  of  the  two-dimensional  case,  the  surface  was  approximated 

by  n  elemental  areas  each  having  constant  surface  current  density  to  develop  a  set  of 

A 

finite  difference  equations.  However,  since  F  (r)  is  composed  of  two  orthogonal 
components  on  each  surface,  the  matrix  equation  will  be  of  order  2n.  The  technique 
was  applied  to  the  problem  of  electromagnetic  scattering  from  a  sphere  due  to  an 
incident  magnetic  field.  The  diameter  of  the  sphere  was  chosen  to  be  .  541 X  (ka  =  1. 7) 
and  72  elemental  areas  described  the  complete  surface.  Because  of  basic  symmetries 
in  the  problem  the  size  of  the  matrix  was  reduced  from  144  x  144  to  36  x  36.  A  more 
complete  description  of  symmetry  conditions  will  be  given  in  subsequent  para^-aphs. 

Figures  12  and  13  are  comparisons  of  the  F£  i0)  and  data  (taken  about  the 
principal  plane)  with  the  claosioal  solutions  of  KL'jg  and  Wu  (Reference  15) . 

The  scattered  fields  are  obtained  from  the  following  equations: 


-■nr  /d "q  fr (St->  |[cog  *  “n  F? 


+  coo 


-  *>FS 
0 


(52) 


H<p  =  hJdsQ~dr  (4“")  |l8in  ®  sUl®  +  cos  9  cos  B  cos  (0  -  *)| 


-  sin  (<p  -  0)  cos  $  F 


9 


(53) 
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The  scattering  cross-sections  are  obtained  from: 


Comparison  of  the  bistatic  cross-section  between  the  "SDT"  and  classical 
results  are  shown  in  Figure  14. 

For  more  complex  geometries  it  is  convenient  to  expand  Equation  8  in  the 
Cartesian  coordinate  system.  Carrying  out  the  indicated  operation  and  resolving  the 
integral  equation  into  sepa.T.te  components  we  obtain, 


2  (K 


n 

oy  p 


Fx  (P)  + 


*4 11 


XX  X 


2  (Ho*  h  •  "oxV  *  *  Fy <p)  +  V=  (lTyx  FX (Q) 


+  TyyFy(«  +  TyIFz(Q)|d.<J 


(55) 


(H  m  -  H  „ 
ox  p  oy 


-  Fz  (P)  + 


Tn 


/JT: 


s 


F  (Q)  +  T  F  (Q)  +  T  F 
zx  xv^  zy  y v^'  zz  z 


(Q)]ds 


H0  (P)  =  H 


OX 


A 

X 


+  H 


oy 


A 

y +  h 


oz 


A 

z 


Tjcx  a  R  (?i  Q)  |(*P  -  Zq>  np  +  (yp  -  yp)  mp 


where 


TXy  =  -  R  (P,  Q)  (Xp  -  Xq)  nip 

^xz  =  "  R  (xp  "  xq)  np 

iy  =-R(P,Q)(yp-yQ)  /p 

Tyy  =  R  (P,  Q)  |  (Xp  -  Xg)  /p  +  (Zp  -  Zq)  np| 

TyZ  =  -  R  (P,  Q)  (yp  -  yp)  np 

Tzx  =  -R<P’<»<zP-2q)  fP 
Tzy  =  'R(p,Q)  (zp‘zq)  mP 

Tzz  “ R  ^p»  ( ^*p  “  xq)  ^p  +  Cyp "  yg)  nip  J 

and 

R  ( P.Q)  -(1  +  Jk|rp  -rQ|) 

If  the  total  aurface  la  approximated  by  n  aubaurfacea  and  we  let  &  Sp  denote  the 
subsurface  at  point  Q,  the  integral  of  Equation  56  nan  be  replaced  by  a  summation 
preparatory  to  numerical  integration. 

2  <Hoy  nP  -  Hoz  mP>  ■  *  F*  (P)  +  Tnt  lTx*  F*  <• 

♦  Txy  Fy  »  +  Txz  Fz  <«!  4S<j 


(56) 


+  Tzy  Fy  '«>  +  Tzz  Fz  <«>!  *SQ  _ 

It  should  be  noted  that  die  summation  applies  to  all  n  surfaces  except  for  P  =  Q 
which  is  the  singularity  term  represented  by  the  single  term  outside  the  summation  on 
the  right  hand  side.  The  left  hand  side  represents  the  incident  wave,  Htac  (P>. 


Each  equation  in  56  represents  n  equations,  (P  *  1  to  n)  and  the  3n 
simultaneous  linear  equations  yield  solutions  of  Fx,  Fy  and  Fz  on  each  subsurface. 

Equation  56  in  matrix  form  is, 

[T||r|  =  |K]  (57) 

where 

J  K  J  *  incident  field  *  2  H4nC  (P)  x  lip 

j  F  J  =  unknown  currents 

(  _  |  Kernel  of  the  integral  equation 

|  1  |  "  (coefficients  of  F„,  F  and  I  In  Equation  56) 

*  y  * 

f  F]  is  readily  obtained  by  multiplying  EqUAtfcm  67  through  by  the  inverse  of  the 
« oefficWnt  matrix, 


-  .3  f 


(58) 


Because  of  the  limited  storage  capacity  of  computers  maximum  accuracy  is 
=ta.ined  by  utilizing  the  symmetry  properties  of  the  scatterer.  For  the  ease  of 
ially  symmetric  bodies  the  number  of  subsurfaces  can  almost  be  quadrupled 
=th  the  same  computer  storage.  When  such  properties  exist  it  is  possible  to  combine 
M  the  coefficients  of  the  symmetrical  terms  and  thereby  reduce  the  number  of  un- 
cvvns .  Care  should  be  taken  to  determine  the  proper  direction  of  the  currents  hence 
— i  signs  of  the  coefficients  before  combining.  For  the  case  of  incident  plane  wave 
zmveiing  along  the  z-axis  and  H  polarized  in  the  y  direction,  the  polarities  of  current 
mponents  on  an  axially  symmetric  body  are  given  in  Table  1  with  the  first  quadrant 
reference  and  the  other  quadrants  given  in  the  conventional  x-y  plane  of  the  right 
sanded  Cartesian  coordinates 

Table  1  Current  Polarities  on  an  Axially  Symmetrical  Body 


Quadrant 

F, 

H 

I 

+ 

+ 

4* 

n 

X 

- 

m 

* 

- 

IV 

l_: _ 

For  the  case  of  plane  symmetry  the  number  of  subsurfaces  could  be  doubled 
applying  the  above  procedures.  With  quadrants  I  and  IV  as  reference,  and  y-z 
!  plane  of  symmetry,  the  polarities  of  the  current  components  are  given  in  Table  2. 

Table  2  Current  Polarities  on  Bodies  with  Plane  of  Symmetry 


Quadrant 
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I&  IV 
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As  an  example  a  spheroid  of  12  x  12  subsurfaces  will  be  considered 


$  =-  o* 
0  =90® 


-  /. 


0  =  0°  A  9  =  90° 


0  =  90° 
9  -  90° 


Sketcn  1  Spheriod  with  12  x  12  Subsurfaces 

In  Sketch  1  the  total  surface  is  divided  into  twelve  sections  in  both  6  and  0 
directions  and  a  set  of  indices  (J,  K)  is  assigned  to  each  subsurface.  The  first  number 
denotes  the  order  in  the  6  direction  and  the  second  number,  the  0  direction.  Table  3 
lists  the  symmetrical  subsurfaces  whose  coefficients  are  combined  to  reduce  the 
computer  storage. 


Table  3  Symmetrical  Subsurfaces  for  Axially  Symmetric  Shapes 


Quadrant 

Symmetrical  Subsurfaces  (K  =  l  to  12) 

I 

U,K) 

(2,  K) 

(3,  K) 

n 

(8,  K) 

(5,  K) 

(4,  K) 

m 

(7,  K) 

(8,  K) 

(9,  K) 

IV 

(12,  K) 

(11,  K) 

(10,  K) 

Once  the  current  on  the  surface  of  a  scatterer  is  determined,  the  fer-fleld 
pattern  can  be  obtained  by  applying  Kirchoff's  integral.  Let  E  ,  E  and  E  denote 

A  y  i* 

the  three  orthogonal  components  of  the  electric  field  at  point  R  *9,  0),  we  have 
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a  v  a  j(kxn,-kyn+k 

EX  =  C  I  Fx^)e  Q  yQ  ^ 


X  Q=1  X 


v  A  i(k  x  +k  y  +k  /.  ; 

EV  =  C  2  Fv^>e  Q  Q  ^ 

•  q«i  y 

E  *  C  2  F,  (Q)  e  W  y  W  ^ 


2  <f?i  z 


where  (x^,  y^,  z^)  =  coordinates  of  the  current  points 

kx  =  k  sin  0  cos  4 
k  =  k  sin  9  sin  4 

y 

k  *  k  cos  9 
z 

C  *  normalization  constant. 

Transforming  into  spherical  coordinates,  we  obtain 

A  A  A  A 

E .  =  Ew  cos  $  cos  4  +  E„  cos  0  sin  #  -  E.  sin  9 
#  »  y  z 

A  A  A 

E^  =  -  Ex  sin  ^  +  Ey  cos  ^ 


Eg  and  E^  as  well  as  Ex,  Ey  and  ig  are  complex  (fliantltles  and  are  expressed  in 
terms  of  real  and  imaginary  components.  When  the  magnitude  and  the  phase  are 
desired  it  is  only  a  matter  of  simple  conversion,  e.g. 


eU 


A‘ 


2*e2 
t  +EI 


,-lEI 
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.4.7W* 


where 


E  =  E  +  i  Et 
R  I 

The  scattering  cross  section  at  r  is  then  given  by 


2.2.4  Application  to  Arbitrary  Three-Dimensional  Geometries 

To  demonstrate  the  versatility  of  the  "SDT"  to  geometries  which  are  not 
tractable  by  classical  procedures,  several  basic  shapes  were  examined.  A  prolate 
spheroid  with  minor  axis  ka  =  1. 7  and  major  axis  kb  =  3. 4,  was  selected  for  com¬ 
parison  with  the  aforementioned  sphere.  In  Figure  15,  propagation  of  the  incident 
magnetic  field  is  along  the  major  axis  and  the  bistatic  cross  section  is  observed 
throughout  the  360  degree  sector.  Considerable  reduction  in  the  monostatic  cross 
section  (in  the  180  degree  direction)  was  achieved.  (Refer  to  Figure  14). 

The  truncated  cone  and  cone-sphere  geometries  are  of  current  interest  since 
they  conform  to  basic  re-entry  vehicle  configurations.  The  bistatic  cross-section 
of  a  truncated  cone  due  to  a  plane  wave  incident  on  the  tip  is  shown  in  Figure  16. 
The  large  back  lobe  is  due  to  the  sharp  discontinuity  at  the  base  of  the  cone.  By 
capping  the  base  with  a  hemisphere  to  elminate  the  discontinuity,  it  1b  possible  to 
suppress  (by  approximately  lOdb)  the  backlobe  contribution  as  shown  In  Figure  17. 

3.0  NORTHROP  NORA  IB '8  "PHYSICAL  OPTICS'’  PROCEDURE 
3.i  General  Fprgwlatipii 

As  indicated  above,  solution  of  the  general  electromagnetic  scattering  and 

diffraction  problem  via  Northrop  Norair's  Source  Distribution  Technique  ("SDT") 
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Is  based  upon  the  explicit  determination  of  the  current  distribution  induced  on  the 
surface  of  the  scatterer  by  the  incident  wave.  The  scattered  fields  anywhere  in  space 
are  in  turn  evaluated  by  integrating  the  current  distribution  over  the  entire  surface. 
The  unique  feature  of  the  "SDT"  is  its  extreme  simplicity  and  flexibility  to  solve 
arbitrary  geometries.  Since  the  scatterer  is  described  in  the  Cartesian  coordinate 
system,  it  can  be  approximated  by  a  finite  number  of  elemental  surfaces  completely 
independent  of  rotational  or  plane  of  symmetry.  Another  advantage  to  choosing  this 
coordinate  system  is  in  the  simplicity  in  evaluation  of  the  scattered  fields  from  the 
current  distribution.  Inasmuch  as  the  three  surface  current  components  on  each 
elemental  surface  are  expressed  in  the  Cartesian  coordinates,  simple  trigonometric 
functions  relate  the  currents  to  the  scattered  fields. 

When  the  characteristic  dimension  of  the  scatterer  lies  in  the  resonance  region 
(.  IX  to  2X)  the  induced  current  distribution  is  obtained  by  numerical  solution  of  an 
integral  equation.  The  integral  equation  is  approximated  by  a  finite  difference  tech¬ 
nique  which  recasts  the  integral  into  a  set  linear  simultaneous  equations  and  is 
solved  by  matrix  inversion  techniques. 

When  the  characteristic  dimension  of  the  scatterer  is  greater  than  2  X,  a 
fundamental  limitation  of  the  matrix  inversion  technique  for  solving  the  induced 
current  distribution  lies  in  the  storage  capacity  and  speed  of  the  digital  computer  and 
therefore  questions  regarding  the  effectiveness  of  the  techniques  to  handle  these 
~*rger  complex  geometries  naturally  arise. 

Indiscriminate  use  of  storage  tapes  to  accommodate  larger  and  larger  matrices 
would  be  seriously  questioned  since  the  rapidly  Increasing  computer  costs  would 
deter  the  analytical  procedures  in  favor  of  empirical  techniques.  Furthermore, 


since  greater  accuracies  are  also  associated  with  larger  matrices  and  hence  higher 
computer  costs,  cost-effectiveness  must  also  be  considered  and  properly  weighed 
with  accuracy  requirements  compatible  with  the  overall  system. 

To  this  end  Northrop  Norair  has  utiliz  e!  the  well-established  procedures  of 
physical  optics  to  determine  the  current  distribution  on  "large"  perfectly  conducting 
8catterers.  The  surface  of  the  scatterer  is  approximated  by  n  elemental  areas 
(expressed  in  Carte  jian  coordinates)  and  the  surface  current  components  are 
evaluated  from  the  tangential  component  of  the  incident  magnetic  field  at  each  sub¬ 
surface. 

A  A  ^ 

J  =  2n  x  H. 
s  t 


where  Jg  =  current  density,  h  *  subsurface  normal, 

/ 

Ht  =  tangential  component  of  incident  magnetic  field  # 
The  scattered  field  is  given  by 


2  IT 

where  k<  j  =  wave  number,  r  distance  from  any  subsurface  to  the  noint 
of  observation  and  die  summation  is  carried  out  over  the  entire  illuminated  area. 


As  can  be  seen  by  the  formulation  above,  a  rigorous  computation  of  KCS  of 
even  the  simplest  shape  is  a  formidable  task,  restricted  to  analytical  surfaces.  The 
results  for  some  of  the  tractable  geometries  are  tabulated  in  Table  4  along  with 
limits  of  applicability  and  proper  formulation  for  regions  of  the  dependent  variables 
(wavelength,  viewing  angle,  etc.).  Generally  the  analyst  must  decompose  the  vehicle 
by  experience  and  good  judgement  to  allow  application  of  catalogued  shapes  and  sum 


I 

to  approximate  RCS.  To  overcome  these  restrictions,  Northrop  Norair  utilized  j 

numerical  procedures  developed  in  part  by  the  previous  "SDT"  program.  The 

surface  of  the  scatterer  is  approximated  by  n  elemental  areas  (expressed  in  Cartesian  | 

coordinates)  and  the  surface  current  components  are  explicitly  evaluated  and  properly 

stored  in  the  computer.  The  scattered  fields  and  RCS  are  then  determined  from  1 

Kirchoff's  integral.  This  computer  program  is  designated  "Physical-Optics"  and  is  j 

applicable  to  surfaces  whose  characteristic  dimensions  are  large  (>  2X). 

3.2  Numerical  Results 

To  test  the  range  of  applicability  of  this  procedure,  the  classical  problem  of 
scattering  from  a  perfectly  conducting  sphere  due  to  an  incident  magnetic  field  was  I 

examined.  Comparison  of  the  backscattering  cross-section  with  classical  results  as 
a  function  of  sphere  circumference  is  shown  in  Figure  18.  The  agreement  is  quite 
good  for  values  of  ka  greater  than  11.  bn  all  cases  the  illuminated  portion  of  the 
sphere  was  approximated  by  2, 820  subsurface  elements.  As  a  matter  of  interest. 

Figure  19  compares  the  bistatic  cross* section  for  a  sphere  (ka  ~  10)  with  the  classical 
series  expansion  results. 

The  flexibility  of  this  "Physical  Optics"  technique  for  the  solution  of  more 
complex  geometries  is  dsmonstr&ted  fay  sxa mining  the  monostatic  scattering  from  a 
large  finite  cylinder.  The  customary  procedure  to  handle  this  geometry  would  be  to 
approximate  it  by  simpler  geometries  for  which  solutions  are  known  (1.  e. ,  the  ends 
by  discs,  cylindrical  body  by  finite  cylinder)  and  then  sum  the  result*  In  an  appro¬ 
priate  manner  to  obtain  the  total  croaa-sectton  of  the  original  shape.  These  result* 
are  included  In  Figures  20  and  21  together  with  measured  data  obtained  from  Northrop 
Norair's  radar  reflectivity  range.  The  discrepancies  in  vicinity  of  45  degrees  are 
due  to  edge  effects  which  are  ignored  in  the  "Physical  Optics"  analysis.  Edges  and 
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surfaces  discontinuities  can  be  handled  by  including  the  "SDT"  procedure  with  "Physical 
Optics"  as  outlined  next. 


4.0  "SDT-PHYSICAL  OPTICS"  PROCEDURES 

The  "Physical  Optics"  approximation  technique  fails  to  describe  properly  the 
currents  on  bodies  in  the  transition  region  (upper  end  of  resonant  region  and  lower 
end  of  physical  optics  region).  Northrop  Norair's  SDT  and  "Physical  Optics" 
techniques  are  suitably  combined  to  determine  explicitly  the  current  distributions 
over  the  entire  surface  and  then  sum  them  to  obtain  the  scatter  fields.  The  combina¬ 
tion  of  the  two  programs  is  identified  as  "SDT -Physical  Optics"  ("SDT-P/O"). 

In  Figure  22,  let  ACB  be  the  illuminated  region  and  ADB,  the  shadow  region. 
The  combination  procedure  simply  applied  the  "Physical  Optics"  techniques  to  a 
portion  of  the  illuminated  area  and  then  determines  the  current  distribution  in  the 
shadow  region  by  applying  die  "SDT"  program,  considering  the  coupling  effects  from 
the  entire  surface. 

The  improvement  of  the  "SDT-Physioal  Optics"  program  over  the  "Physical 
Optics"  technique  for  bodies  in  the  transition  region  is  clearly  demonstrated  in 
Figure  23.  It  is  clear  that  a  more  aceurate  description  of  the  currents  in  the  shadow 
region  has  been  achieved  in  the  combination  process  resulting  in  more  accurate 
description  of  the  scattered  fields. 

5.0  EXPERIMENTAL  APPARATUS 

An  image-plane  technique  was  utilized  to  measure  the  current  distribution  and 
scattered  fields.  The  general  instrumentation  followed  along  the  lines  of  Kodis 
and  Wetzel  (Reference  18)  who  successfully  applied  it  to  two-dimensional  surfaces. 


A  few  modifications  to  their  procedure  was  necessary  to  accommodate  the  three- 
dimensional  geometries  described  in  this  study. 

An  image-plane  was  constructed  of  1/16  inch  aluminum  sheeting  and  electrically 
bonded  to  form  a  continuous  16-foot  x  24-foot  ground  plane.  Flatness  was  held  to 
1/8  inch  over  the  entire  surface.  Located  in  the  center  of  the  ground  plane  was  a 
circular  opening  15  inches  in  diameter  The  well  beneath  the  cut-out  housed  the 
measuring  and  positioning  apparatus.  An  aluminum  plate,  together  with  the  position¬ 
ing  device,  was  fitted  into  the  hole  flush  with  the  surface  and  secured  from  below. 
Hemi-surfaces  of  various  shapes  investigated  were  fabricated  and  fastened  to  the 
aluminum  plate  (See  Figure  24). 

For  the  current  distribution  measurements,  a  magnetic  probe  mounted  on  a 
copper  band  was  carried  over  the  surface  for  various  elevation  angles.  This, 
together  with  the  option  of  rotating  the  circular  plate  in  azimuth  afforded  the  system 
with  two  degrees  of  freedom.  Amplitude  and  phase  meausrements  for  both  elevation 
and  azimuth  positions  were  thus  obtained. 

Plane  wave  excitation  was  approximated  by  an  H-plane  flared  horn  fastened 
at  one  end  of  the  image-plane.  Since  the  mode  of  propagation  over  the  image-plane 
was  TM,  all  measurements  were  conducted  for  only  vertical  polarization. 

Calibration  of  the  system  was  performed  on  conducting  spheres  for  which 
formal  solutions  are  well-known  and  extensively  tabulated.  Results  of  these  measure¬ 
ments  are  shown  in  Figures  13  and  14.  The  primary  source  of  error  was  from  re¬ 
flecting  edges  of  the  ground  plane.  To  minimize  the  edge  effects,  absorbing  material 
was  laid  and  secured  along  the  entire  periphery.  (See  Figure  25).  Bistatic  reflectivity 
measurements  were  also  conducted  on  the  image-plane  utilizing  an  interferometer 
technique.  In  operation  the  bridge  is  initially  nulled  without  the  test  specimen,  then 
with  it  in  place;  the  recorded  difference  is  the  measure  of  the  scattering  cross* section. 


The  Image  plane  is  also  equipped  with  a  rotator  and  transmitter  mount  to 
function  as  a  ground  plane  reflectivity  range  as  Illustrated  In  Figure  26.  A  more 
complete  description  of  the  measurement  apparatus  and  test  procedures  will  be 
published  in  a  forthcoming  report . 

6.0  CONCLUSIONS 

The  results  of  this  study  indicate  that  a  useful  and  versatile  technique  has  been 
developed  for  the  solution  of  electromagnetic  scattering  from  arbitrary  shapes.  It  is 
capable  of  solving  two-  and  three-dimensional  scattering  problems  for  which  classical 
solutions  are  not  available.  For  those  simple  geometries  for  which  exact  solutions 
are  available  for  comparison,  excellent  agreement  was  found  in  all  cases.  For  the 
other  shapes  the  "SOT''  solutions  correlate  well  with  the  measured  data. 

Since  the  scatterer  is  described  in  the  Cartesian  coordinate  system,  it  can  be 
approximated  by  a  finite  number  of  elemental  surfaces  completely  independent  of 
rotational  or  plane  of  symmetry.  Another  advantage  of  choosing  this  coordinate 
system  is  the  simplicity  In  the  eval  uation  of  the  scattered  fields  from  the  current 
distribution.  As  the  three  surface  current  components  on  each  elemental  surface 
are  expressed  in  the  Cartesian  coordinates,  simple  trigonometric  functions  relate 
the  currents  to  the  scattered  fields  hi  their  appropriate  designations. 

Although  developed  primary  for  the  "resonant"  frequency  range,  the  "SDT" 
program  has  been  suitably  combined  with  large  body  approximation  procedures, 

"SDT- Physical  Optics,"  handle  larger  geometries  with  greater  ease  and  flexibility 
then  present-day  approximation  procedures. 
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FIGURE  1  CURRENT  DISTRIBUTION  ON  CIRCULAR  CYLINDER,  Ht,ne  (VERT.  POL.)  (ko  ■  1.6) 
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FIGURE  2  CURRENT  DISTRIBUTION  ON  CIRCULAR  CYLINDER,  Exi,,e  fHOR.  POL.)  (ka  ■  1.6) 
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FIGURE  3  SCATTERING  CROSS-SECTION  OF  CIRCULAR  CYLINDER,  (VERT.  POL.)  (ko  -  1.6/ 
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FIGURE  4  SCATTERING  CROSS-SECTION  OF  CIRCULAR  CYLINDER,  E,1"*  (HOR.  POL.)  (ko  *  1.6) 
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IN  DEGASES  AMPLITUDE 


CURRENT  AMPLITUDE  DISTRIBUTION  ON  ELLIPTIC  CYLINDER,  HIinc(VERT.  POL.) 
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FIGURE  6  CURRENT  PHASE  DISTRIBUTION  ON  ELLIPTIC  CYLINDER,  H,1"*  (VERT.  POL.) 
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FIGURE  7  SCATTERING  CROSS-SECTION  OF  ELLIPTIC  CYLINDER,  Hfine  (VERT.  POL.) 
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FIGURE  8  AIRFOIL  SECTION  NO.  1 
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AMPLITUDE  AMPLITUDE 


SOURCE  DISTRIBUTION  TECMN'QUE  (“SOT”) 
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1  MEASUREMENT 


o  I  - - a - - - - ;  i - 1 . . . . . . .  ,i 

C  30  M  «0  120  ISO  1*0 


FIGURE  12  CURRENT  DISTRIBUTION  ON  SPHERE  F^$)  COMPONENT,  (ko  *  1.7) 


18  CURRENT  DISTRIBUTION  ON  SPHERE  Fjt)  COMP"*': NT,  (ko  *  1.7) 
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FIGURE 


FIGURE  19  BISTATIC  CROSS-SECTION  FOR  PERFECTLY 
CONDUCTING  SPHERE  (ko  - 10)  H-PLANE 


20  MONOITATIC  CROSS'SCCTION  OF  CIRCULAR  CYLINDER 
(R*  U1\,  L-2.42M 
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"PHYSICAL  OPTICS" 
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FIGURE  23.  "SDT  -  PHYSICAL  OPTICS"  BBTATIC  CROSS  SECTION  OF  SPHERE  (ks  *  4 
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FKJURE28.  GROUND  PLANE  REFLECTIVITY  RANGE 
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TABLE  1  RADAR  CROSS  SECTION  ( r  )  OF  SOME  SIMPLE  SHAPES 


_  _  _ _ fmmmr  mum 

TABLE  1  RADAH  CROSS  SECTION  (a)  OF  SOME  SIMPLE  SHAPES 

(CONTINUED) 


SHAPE 

PERTINENT 
DIMENSIONS  (AND 
POLARIZATION  >• 

viewing 
ORIENTATION 
i  AND  BODY 

WAVE  LENGTH 
CONDITIONS)** 

RADAR  CROSS  SECTION 

0 

Paraboloid 

EBB 

At  angle  «.  from  Axis 

2  4 

4tt  p  sec 

Wire 

Euler's  constant  =  >  * 
i.7* 

Radius  of  wire  -  a 

Length  of  wire  -  L 
Polariiatton  at  angle  # 
from  parallel  to  axis 

At  angK  .  to  Axis 

Of 

tin  L  \ 
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